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Resonant x-ray scattering (RXS) experiment has been performed for the (3 1.5) superlattice 
reflection in the antiferroquadrupolar and antiferromagnetic phase of DyB2C2- Azimuthal-angle 
dependence of the resonance enhanced intensities for both dipolar (El) and quadrupolar (E2) 
resonant processes has been measured precisely with polarization analysis. Every scattering channel 
exhibits distinctive azimuthal dependence, differently from the symmetric reflection at (0 0.5) 
which was studied previously. We have analyzed the results using a theory developed by Lovesey et 
al., which directly connects atomic tensors with the cross-section of RXS. The fitting results indicate 
that the azimuthal dependences can be explained well by the atomic tensors up to rank 2. Rank 
3 and rank 4 tensors are reflected in the data very little. In addition, The coupling scheme among 
the 4/ quadrupolar moment, 5d ortitals, and the lattice has been determined from the interference 
among the Thomson scattering from the lattice distortion and the resonant scatterings of El and 
E2 processes. It has also been established from the RXS of the (3 1.5) reflection that the canting 
of the 4/ quadrupolar moments exists up to Tq. We also discuss a possible wavefunction of the 
ground state from the point-charge model calculation. 

PACS numbers: 75.25.+Z, 61.10.Eq, 71.20.Eh, 75.40.Cx 



I. INTRODUCTION 

Anisotropic charge distributions of localized 4/ elec- 
trons play important roles in magnetic properties of rare- 
earth compounds through quadrupolar interactions and 
orderings since the magnetic moments are strongly cou- 
pled with the quadrupolar moments by the spin-orbit 
interaction. Unusual magnetic phase diagrams and field 
induced antiferromagnetic structures typically observed 
in CeBg, TmTe, DyB2C2, etc., cannot be understood 
without considering the quadrupolar moments i^i&ii^ In 
CeBg even an octupolar moment is necessary to explain 
the experimental observations £L£ What is more exotic, 
in NpC>2, it is claimed that the IVtype octupolar or- 
der is realized without a magnetic dipolar moment £ It 
is therefore very important to study various behavior of 
quadrupolar and higher rank multipolar moments. 

Recently, a series of tetragonal rare-earth compound 
of RB2C2 (R=Tb, Dy, Ho) has been studied extensively 
as a typical system which exhibits antiferroquadrupo- 
lar (AFQ) orderings and complex magnetic struc- 
tures ^iSaii^iSii 3 . Among the three compounds DyB2C2 is 
the most typical AFQ system with the highest Tq at 24.7 
K and a well separated T/v at 15.3 K, indicating that the 
quadrupolar interaction is stronger than the magnetic in- 
teraction. In H0B2C2 and TbB2C2 the two interactions 
are comparable in magnitude and compete with each 
other, leading to complex magnetic structures and phase 
diagrams^ The magnetic structure of DyB2C2 in the 
antiferromagnetic (AFM) phase below T/v determined by 
Yamauchi et al. is shown in Fig.^ 5 The moments along 
the c axis are almost perpendicular with each other and 
those in the ab plane are antiparallel but canted from 



the (1 1 0) directions. This structure is characterized by 
four fc-vectors: fe 1= (l 0), fc 2 =(l 0.5), fc 3 =(0 0), 
and fc4=(0 0.5). This unusual magnetic structure can 
be understood naturally by considering the underlying 
AFQ order. 

Observation of the AFQ order in DyB2C2 has been 
performed to date by resonant x-ray scattering (RXS). 
By tuning the energy of the incident photon to an ab- 
sorption edge of the element in study, a core electron 
is promoted to an unoccupied valence shell, forming an 
intermediate state. Secondary photon is emitted when 
the electron returns to the core state, and the scattered 
photon contains information of the valence shell. At the 
Liu absorption edge of the rare-earth elements the in- 
termediate state is described by a 2p 3 / 2 core hole and 
an additional electron in the 5d shell (El process) or 4/ 
shell (E2 process). 

After the first observations of RXS in DyB2C2 by Hi- 
rota et al. and Tanaka et aZ. ) 14 i 15 we made a detailed 
analysis on the observed reflections using a theory de- 
veloped by BlumeiiSiii The experimental results are well 
explained by the model of the AFQ structure shown in 
Fig. n In Blume's theory, the scattering-factor tensor is 
devided into symmetric and asymmetric parts, which are 
interpreted in Ref. as quadrupolar and magnetic ori- 
gin, respectively. The scattering amplitudes calculated 
for each fc-vector, polarization channel, and transition 
process (El and E2), are mostly consistent with the ob- 
served results. 

In Ref. 0, we suggested that the canting angle of the 
quadrupolar moment above Tjv could be zero; i.e., they 
align along the [1 1 0] directions. This statement came 
from the experimental result that the (1 2.5) resonant 




FIG. 1: Top view of the crystal structure of DVB2C2 
(PA/rnbrn, a = 5.341 A, c = 3.547 A at 30 K). Arrows in- 
dicate the magnetic moments in the AFM phase. Ellipses 
indicate the schematic views of pancake-like charge distribu- 
tions of Dy below Tq . Their flat planes are perpendicular to 
the magnetic moments. The B-C layers are shifted from the 
Dy layers by z — 0.5. 

reflection for the a- n' channel seemed to disappear above 
Tat. The calculated scattering amplitude of this reflec- 
tion in the AFQ phase contains the symmetric part of 
the scattering-factor tensor multiplied by sin 2a, where 
a is the canting angle from the [1 1 0] directions. Since 
the former factor is definitely not zero in the AFQ phase, 
the latter factor a must be zero if the experimental re- 
sult is correct. This is an important problem when we 
consider the origin of the canting in detail; one is an 
inter-ionic quadrupolar interaction mediated by the con- 
duction electrons and the other is a quadrupole-strain 
interaction by the tilts and shifts of the surrounding B- 
C octagon. Note that the B-C octagons around Dy- 
1 and Dy-2 are oppositely tilted about the c axis. If 
a = 0° above Tjy, we would come to a statement that 
the inter-ionic quadrupolar interaction is much stronger 
than the quadrupole-strain interaction and prefer 90° ar- 
rangement along the c axis and parallel arrangement in 
the ab plane. 

Lattice distortion below Tq characterized by the prop- 
agation vector of (1 0.5) is also important in DVB2C2. 
Although this lattice distortion should be compatible 
with the AFQ order, we tried to discuss a possibility of 
the motion of Dy atoms in the previous papersii^i^ We 
compared the structure factors calculated for two models 
of the subgroup of P4/mbm; in one model the Dy atoms 
are shifted (PA/mnc, No. 128) and in the other the B-C 
atoms are shifted (P^/mnm, No. 136), both along the 
c axis. In P^/mnm B-l and B-3 in Fig. ^ shift upward 
and B-2,4 shift downward. Motions of C atoms have two 
possibilities: C-1, 3 up and C-2,4 down or C-2,4 up and 
C-1, 3 down. We calculated the former case in Ref. ITfjl 
and claimed that the observed intensity of the (1 2.5) 
reflection cannot be obtained. However, the latter case 
explains the observed intensity by assuming the shift of 
about 10 _3 c. Later on, Adachi et al. performed non- 
resonant x-ray diffraction experiment and observed (0 1 
half-integer), (11 half- integer), and (1 2 half- integer) re- 
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FIG. 2: Left: Configuration of the incident and scattered x 
rays in the laboratory coordinates XYZ. Right: Configura- 
tion of the reciprocal lattice space in the laboratory coordi- 
nates at the azimuthal angle \& = 0°. The sample rotation is 
counter-clockwise about Z with increasing ^> . 



flections. The intensities are well explained by assuming 
the B-C motion of the latter case^ The motion of Dy in 
PA/mnc forbids the (11 half- integer) reflections, which 
contradicts the experimental result. It should be noted 
that the validity of the analysis of the RXS data is not 
affected because it assumes the site symmetry of 2/m 
which is compatible with the space group P^/mnm. 

In the previous works of RXS in DyB2C2, the (0 half- 
integer) reflections for the a-a' channel in the E2 process 
were hidden in the tail of the El resonance, while those 
for the cr-7r' channel were well resolved i£il£ilfL Recently, 
Tanaka et al. made a systematic RXS study of the (0 
half-integer) reflections and analyzed the data using the 
structure factor of atomic tensors calculated by Lovesey 
et alm^^ They ascribed the disappearance of the E2 res- 
onance in the a-a' channel to the idea that the rank 2 
(quadrupole) and rank 4 (hexadecapole) tensors cancel 
out in the structure factor. Although they claim this as 
the observation of the rank 4 tensor, it should be checked 
by more detailed experiments because such a high rank 
multipole should be difficult to observe experimentally. 
We have shown in another paper that the rank 4 multi- 
pole is not necessary if we consider the interference be- 
tween El and E2 resonancesiSi Further evidence that 
support this statement will be presented in this paper. 

In the present study, we have performed RXS experi- 
ment on the (3 1.5) reflection to investigate the above 
problems; one is the canting angle in the AFQ phase and 
the other is the possibility of observing the higher rank 
atomic tensors. We have analyzed the data using a the- 
ory developed by Lovesey et at, which directly connects 
atomic tensors of different ranks with the cross-section 
of RXSi 22 i 23 i 24 i 25 A theory by Blume^ or Hill and Mc- 
Morrow, 2 - basically gives the same result as the one by 
Lovesey et al. However, physical meanings of the param- 
eters that appear in the formulae are not as clear as those 
in Lovesey's formulae. 
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II. EXPERIMENT 

A single crystal of DyB2C2 was grown by the Czochral- 
ski method with a tetra-arc furnace. A flat (2 1) surface 
with its area about 1x2 mm 2 was obtained by rub- 
bing the crystal with sandpaper and then polishing the 
surface with fine emery paper. RXS experiments were 
performed on a four-circle diffractometer at the beam- 
line 16A2 of the Photon Factory, KEK. The incident en- 
ergy was tuned near the Dy Lni edge. The sample was 
mounted in a closed-cycle helium refrigerator so that the 
(2 1) plane was normal to the </> ax is of the diffrac- 
tometer. The configuration of the sample and x rays are 
illustrated in Fig. [21 The azimuthal-angle scan was per- 
formed by rotating the <j> axis of the diffractometer. The 
mosaic width of the (2 1) fundamental reflection was 
about 0.06°. 

The incident x ray is almost linearly polarized with 
its electric field perpendicular to the scattering plane (<r 
polarization). To separate the a' (unrotated) and tt' 
(90° rotated) components of the diffracted beam, PG 
(0 6) reflection was used as an analyzer, for which 
the scattering angle becomes 90.7° and almost perfect 
analysis (cos 2 90.7° = 1.5 x 10~ 4 ) is achieved. However, 
the contaminations caused by the unpolarized compo- 
nent of the incident beam cannot be neglected. The to- 
tal contamination, estimated from the intensity ratio of 
^rotated /-f unrotated of the fundamental (2 1) reflection, 
was about 0.3%. 



III. RESULTS AND ANALYSIS 
A. Energy dependence 

Figure shows the energy dependences of the (3 1.5) 
reflection at T = 12 K in the AFM phase and at T = 20 
K in the AFQ phase. The data have been corrected for 
the absorption coefficient that was determined from the 
fluorescence spectrumiSi The disagreement of the inten- 
sities at energies well below and above the edge is due 
to the imperfect correction. Large nonresonant contribu- 
tion of the Thomson scattering from the lattice distortion 
is observed. The data are normalized to the intensity at 
7.72 keV and ¥ = 180° because this Thomson scattering 
does not exhibit azimuthal-angle dependence, which has 
been checked by comparing with the azimuthal scan of 
the (2 1) fundamental reflection. 

An azimuthal-angle dependence is clearly observed 
around the edge. It is remarkable that a large resonant 
enhancement is observed at 7.782 keV, which is ascribed 
to the E2 resonance. The enhancement of the El res- 
onance at 7.792 keV is not clearly seen because the en- 
ergy dependence of the nonresonant component due to 
the absorption is much stronger than that of the El res- 
onance. However, the observed azimuthal dependence 
around 7.792 keV may be ascribed to the resonance of 
the El process. 
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FIG. 3: Incident energy dependences of the intensity of the 
(3 1.5) reflection for the a-a' channel at $ = 180° and 
270°. The energy dependence of the absorption coefficient 
determined from the fluorescence spectrum is shown in the 
top figure. 
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FIG. 4: Incident energy dependences of the intensity of the (3 
1.5) reflection for the a-n' channel at * = 180° and 270°. 



Figure 0] shows the energy dependences for the (j-tt' 
channel. Contamination of the unrotated scattering has 
been subtracted and the data are normalized to the non- 
resonant intensity at 7.72 keV and ¥ = 180° for the a-a' 
channel. At 20 K, only the El resonance is observed at 
¥ = 180° and no signal is observed at ¥ = 270°. At 10 
K, the spectrum at ¥ = 180° exhibits large enhancement 
with a long tail to the lower energy side because of the in- 
terference between the nonresonant magnetic scattering 
and the resonant scatterings of El and E2 processes^ 
The spectrum at ¥ = 270° exhibits only the E2 reso- 
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FIG. 5: Azimuthal-angle dependences of the integrated inten- 
sity of the (3 1.5) reflection for the El process. Solid lines 
are the fits described in the text. Dotted and dashed lines in 
the a-n' channel indicate the rank 1 and rank 2 contributions, 
respectively, to the fitting result at 12 K. 

nance. In both a-a' and cr-7r' channels, the resonant fea- 
tures clearly exist still at 20 K in the AFQ phase, which 
indicate that the the quadrupolar moment is still canted. 
The result is in accord with the site symmetry of 2/m 
and a symmetry analysis recently performed by Zaharko 
etalM 



B. Azimuthal-angle dependence 

Figure [3] shows the azimuthal-angle dependence of the 
integrated intensity for the El process measured at 7.792 
keV. The data are normalized to the nonresonant inten- 
sity at 7.72 keV for a-a' . In spite of the experimental con- 
dition that 7.792 keV corresponds to the minimum posi- 
tion of the energy spectrum in Fig. |31 because of the ab- 
sorption, the intensity of the a-a' channel varies with the 
azimuthal angle with 180° periodicity. It is noted that 
the amplitude of the azimuthal dependence for the a-a' 
channel is much larger than that for the a-n' channel. 
This is due to the interference with the strong Thomson 
scattering in the a-a' channel, which will be explained in 
the next subsection. 

Figure shows the azimuthal-angle dependence for 
the E2 process measured at 7.782 keV. The data are 
again normalized to the nonresonant intensity at 7.72 
keV for a-a'. Apparently, the azimuthal dependence is 
completely different from that of the El process. This 
fact also supports the ascription of this resonance to the 
E2 process, which measures the 4/ state directly. This 
is different from the case in YTiC>3 where the azimuthal- 



FIG. 6: Azimuthal-angle dependences of the integrated inten- 
sity of the (3 1.5) reflection for the E2 process. Solid lines 
are the fits described in the text. Dotted and dashed lines in 
the a-n' channel indicate the rank 1 plus rank 3 (magnetic) 
and rank 2 (quadrupolar) contributions, resspectively, to the 
fitting result at 12 K. The dotted line in the a-a' channel 
indicate a fit including the rank 4 contribution. 

angle dependences of the RXS intensity at the pre-edge 
and the main-edge are the same and the pre-edge reso- 
nance is ascribed to the El process^ 

C. Analysis 

1. Azimuthal-angle dependence 

In order to analyze the energy and azimuthal-angle 
dependences, we have calculated the structure factors of 
atomic tensors following a theory developed by Lovesey 
et al22£££&£^ The structure shown in Fig. ^ is assumed. 
The canting angle (f> of the magnetic moments, which 
is assumed to coincide with the principal axis of the 
quadrupolar moment in the AFQ phase, is fixed at 73° 
as determined in Ref. 0. Note that </> is an angle from the 
[1 0] axis in this paper. The spherical tensor (Tq) is 
defined in the local ionic coordinates xyz as indicated in 
Fig. 2] The direction of the magnetic moment is taken 
as the local x axis. 

Since the Dy ion is in the local symmetry of 2/m in 

the AFQ phase, {Tq K ^) with K =even has only the com- 
ponents of q = ±2 and the relation (Tg ) = (T^P) is 
satisfied. Both (T 2 (2) ) and (T 2 (4) ) are real, and (T 2 (2) ) rep- 
resents the ^(x 2 — y 2 ) type quadrupolar moment and 

(T 2 (4) ) represents the ^{x 4 - y 4 - f (x 2 - y 2 )r 2 } type 
hexadecapolar moment. With respect to the odd rank 
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FIG. 7: Calculated azimuthal-angle dependences of the ten- 
sor structure-factors of the (3 1.5) reflection for the El 
resonance. The dashed line for the rank 1 tensor indicates 
the structure factor of the other domain. The domain only 
reverses the sign for the rank 2 tensor. The squares of these 
factors are compared with the data in Fig. 03 
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FIG. 8: Calculated azimuthal-angle dependences of the ten- 
sor structure-factors of the (3 1.5) reflection for the E2 
resonance. The dashed lines for the odd rank tensors indicate 
the structure factors of the other domain. The domain only 
reverses the sign for the even rank tensors. The squares of 
these factors are compared with the data in Fig. |S| 



tensors, since we assume the dipolar moments lie in the 
ab plane, only (Tq K ^) with q =odd is allowed and the 
relation (T^ K) ) = -(T [K q ] ) is satisfied. Then, for K = 3, 

only |(x 3 - fxr 2 ) = Tf u type and ^x{y 2 - z 2 ) = T$ u 
type octupolar moments are allowed. 

The structure factors of the atomic tensors for the El 



and E2 processes are obtained by summing and 
^( S2 ) f or f our Dy j orls i n the unit cell together with 
the phase factors as described in the Appendix. The 
azimuthal-angle dependences of the calculated structure 
factors F^ 1 ) and F^ E2 "> are shown in Figs. and |SJ re- 
spectively. The constant factors A and B are not consid- 
ered. It should be noted that (T { q K) ) in F iE ^ and F iE ^ 
represent atomic tensors of the 5d and 4/ shells, respec- 
tively, and it is written as (Tq K ^) 5 d or (T ? (if) ) 4/ when we 
specify the shell. 

For even rank tensors with K =2 and 4, the coeffi- 
cients of (T^) and (T2 ), respectively, are shown. For 
odd rank tensors with K =1 and 3, the coefficients of 
i(J x ) and i(Tf n ), respectively, are shown. The squares 
of these structure factors are compared with the experi- 
mental results. The fitting results are shown by the solid 
lines in Figs. and At 20 K in the AFQ phase only 
the even rank tensors contribute to the scattering and at 
10 K in the AFM phase the odd rank tensors of magnetic 
origin are superimposed. Since the odd rank tensors ex- 
hibit the domain structure, the intensities from the two 
domains are averaged in the fitting procedure. 

With respect to the El process, the data at 20 K are 
fitted with the structure factor of rank 2 only. The data 
for both polarization channels are well reproduced by the 
calculation, although there is a slight difference in the a- 
a' channel around $ = 180°. At 12 K, the structure 
factor of rank 1 is added in the fitting, and the data 
for <7-7r' channel are well reproduced. The data for a- 
a' channel still exhibit the azimuthal dependence of the 
rank 2 tensor because there is no contribution from the 
rank 1 tensor. 

With respect to the E2 process, the data at 20 K are 
well reproduced by the structure factor of rank 2 only. 
Although the rank 4 tensor should in principle contribute 
to the scattering, the data shows only the feature of the 
rank 2 tensor. The present scattering vector of (3 1.5) 
has different azimuthal-angle dependences for rank 2 and 
rank 4 tensors. Then, the introduction of the rank 4 
tensor only disturbes the fitting results in Fig-EJ For a-a' 
the fitting seems to be slightly improved by introducing 
the rank 4 tensor about 15% as large as the rank 2 tensor. 
The result is shown by the dotted line, but this is within 
the experimental error. It is remarkable that the data 
for er-7r' at 20 K are almost perfectly reproduced only by 
the rank 2 tensor. 

The data at 12 K for the E2 process involve even and 
odd rank tensors. However, for the a-a' channel, the 
data were fitted only by the rank 2 tensor because it was 
not possible to discriminate the azimuthal dependence 
of the odd rank tensors. The data for the cr-7r' channel 
were fitted roughly by the rank 2, rank 1 and rank 3 
tensors. However, the intensity at 7.782 keV for the a-ir' 
channel involves relatively large contribution from the in- 
terference between the nonresonant magnetic scattering 
and the El resonance. Then, the data strongly reflects 
the azimuthal dependence of the El resonance, and it is 
difficult to discriminate the pure E2 resonance. 
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FIG. 9: Simulation of the energy and azimuthal-angle depen- 
dences of the (3 1.5) reflection near the absorption edge, con- 
sidering the interference among the Thomson scattering and 
the resonances of El and F2 processes. The lattice does not 
contribute to the a-n' channel. The parameters are taken to 
be F lat = 20, F< E2 > = F c ( a t 2) x (-0.1), F< B1 > = x 0.033, 

r B i = V E 2 = 6 eV, A B 2 = 7.78 keV, and A B i = 7.79 keV. 



2. Interference effect 

We have so far analyzed the azimuthal-angle depen- 
dences of the El and E2 processes independently. With 
respect to the data at 20 K in the AFQ phase, only the 
rank 2 tensor is enough to explain the data. However, we 
should be careful about the interference between El and 
E2 resonances. It is not obvious if the intensity at 7.782 
and 7.792 keV reflects only the E2 and El resonance, 
respectively. In the (0 half-integer) reflection the inter- 
ference between El and E2 resonances affects the energy 
spectra for a-a' and a-ir' channels as we have shown in 
a former paper^i Furthermore, in the a-a' channel the 
nonresonant Thomson scattering from the lattice distor- 
tion also interfares. We have to consider 
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(1) 



for the <j-g' channel. F\ at is not necessary for the a-ir' 
channel. r\ and T2 are the mixing parameters to model 
more realistic oscillators of El and E2 resonances. In the 
present analysis, however, they are both fixed at unity. 
This point will be discussed in the next section. 

Figure [5] shows a simulation of the energy and 
azimuthal-angle dependences in which the interference 
effect is considered. The calculated structure factors in 
Figs. □ and il are used. The features at 20 K in Figs. [3 El 
03 and El are well explained, although the El resonance 
Fig. |21 at ^ = 270° is hidden in the strong dip due to 
the absorption effect. With respect to a-a' , the scale of 
the longitudinal axis is much larger than that of a-ir' be- 
cause the Thomson scattering from the lattice distortion 
interferes. 

In order to reproduce the experimental results, the 



signs of (T 2 

.(2) 



(2) 

4/ and (T 2 )§d must be opposite and the 



signs of (T^^/ and Fi at must also be opposite. In Fig.0 



we assumed (T 2 (2) ) 4/ = -0.1 and (T^') 5d = 0.033. The 
opposite sign between the two tensors is consistent with 
the result of the RXS of the (0 2.5) reflection in Ref.|2l| 
and the relation with i*i at is newly obtained in the present 

( 2) 

reflection of (3 1.5). We note that the values of (T 2 } 4 / 

( 2) 

and (T 2 '} 5 d used here are slightly different from those 
obtained in Ref. EH for the (0 2.5) reflection. This 
is because the energy dependence of the atomic scatter- 
ing factor is not common, but depends on the reflection 
points, as observed, e.g., in YTiOaiSi 



,(2) 



The matrix element of O 



(2) _ V3, 



■J*) 



in the local 



—if / 2 

'2 — 2 V x 

coordinate becomes positive for the charge distribution 
of the 4/ electrons shown in Fig. ^ which is extended 
along the y direction. However, since the reduced matrix 
element (r.m.e) of for the 4/ electrons of Dy 3+ is 

(2) 

negative, (T 2 ) 4 / becomes negative. Concerning the 5d 
orbital, we should consider yz and zx orbitals which are 
degenerate at T > Tq . Below Tq the degeneracy is lifted 
by the quadrupolar moment of the 4/ orbital through the 
d-f Coulomb interaction or by the lattice distortion of B 
and C atoms through the crystal field effect. Both have 
the local symmetry of 2/m, leading to the appearance of 
K ">■ "■■ 7 • The above analysis indicates that the 

h(2)\ 



(T^') 5d below T c 



sign of (T 2 ') 5 d is positive. Calculation of (T 2 ') 5 d using 
a formula in Ref. |2J shows that the zx orbital gives pos- 
itive and yz orbital gives negative values with the same 
magnitude. This indicates that the zx component in the 
occupied state increases in the AFQ phase. Therefore, 
the charge distributions of 4/ and 5d orbitals are orthog- 
onal with each other. 

Finally, the requirement that the sign of F\ at must be 
positive indicates that the B-C parallelogram of B-2, B-4, 
C-l, and C-3, is shifted up and that of B-l, B-3, C-2, and 
C-4, is shifted down in the layer of z = 0.5; the direction 
of the shifts is opposite in the layer of z = 1.5. In the 
opposite case the sign of Flat becomes negative, which is 
inconsistent with the experimental result. Then, the B-C 
atoms along the direction where the 4/ charge distribu- 
tion is extended move away from the Dy ion, while those 
along the direction where the 5d charge distribution is 
extended move toward the Dy ion. 



IV. DISCUSSIONS 

A. multipolar moments 

The experimental results exhibit clear evidence of the 
rank 2 tensor, but the rank 4 tensor is little reflected in 
the data. Let us discuss the 4/ wavefunction in the AFQ 
phase and the multipolar moments. In order to simulate 
the AFQ phase, we calculated a set of eigenfunctions of 
a crystal field Hamiltonian by the point charge model. 
Since the entropy at Tq is i?ln4^ the ground state must 
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rankO rank 0+2 rank 0+2+4 rank 0+2+4+6 




FIG. 10: Calculated charge distribution of the 4/ electrons 
of a Dy ion in the AFQ phase. The figure is shown in the 
ionic coordinates and the x-axis corresponds to the direction 
of the magnetic moment in the AFM phase. Asphericity is 
emphasized by subtracting the spherical charge distribution 
of a Gd ion with seven 4/ electrons. 



be a pseudo-quartet where the first excited doublet has 
much smaller energy than the second one. In addition, 
since the dipolar moment has large components within 
the ab plane, the main components of the eigenfunctions 
of the pseudo-quartet must have small J z . The eigen- 
functions which satisfy these conditions were obtained by 
appropriately selecting the effective point charges. Next, 

(2) 

in the local coordinates, an 2 -type perturbation was 
added to make the eigenstates in the local symmetry of 
2/m. Linear combinations within the pseudo-quartet are 

(2) 

made and the eigenfunctions which diagonalizc 2 were 
obtained. The angular part of the charge distribution cal- 
culated for thus obtained ground doublet is demonstrated 
in Fig. ED The charge distribution is expressed by the 
even rank spherical harmonics up to rank 6, and the fig- 
ure is shown with increasing the rank step by step. This 
eigenfunction has a magnetic dipolar moment of 9.5/is 
along the x direction, which is consistent with the exper- 
imental result i^i 

The higher rank multipolar moments 2 with even 
K are simultaneously diagonalized. The expectation val- 
ues of the operator equivalent 2 K ^ for K = 2 and 4 are 
listed in Table U together with the r.m.e for the atomic 
tensor T^ K > that appear in the scattering factor of the 
E2 resonance^ The expectation value of (T^^if is ob- 
tained by multiplying the r.m.e and (0 2 K ^). The Stevens 
factors dfj^ are also listed^ The absolute values of the 

multipolar moments are equal to 6j ■ (0 2 •*)• Note that 
the value for rank is nine, which is the total number of 
the 4/ electrons. 

This calculation shows that the rank 4 moment (= 
—0.071) is not negligible in comparison with the rank 2 
moment (= —0.225). This can be observed in Fig. 1101 in 
the difference of the charge distribution including up to 
rank 2 and rank 4 moments. An estimation of {T 2 K ^)$f 
using the r.m.e in Table [3 gives much larger value of 

(T 2 (4) ) 4/ (= 242) than (T 2 (2) ) 4/ (= -18.5). However, 
the present experimental results, when compared with 
the calculated structure factors in Fig. |H| indicate that 
(T^ 4 ^/ is much smaller than (T, 2)4/. This is inconsis- 
tent with the above simulation. 



TABLE I: Matrix element of the operator equivalent Oi, , 
reduced matrix element of T^ K ' for the E2 resonance at the 
Lin edge (j = |), Stevens factor, and the multipolar moment 
of the calculated wavefunction in Fig. 1101 for rank 2 and 4. 
The definition of Of ] follows Ref. H 



K 


(Of ) 


(J\\T(K:j)\\J) 




eT(o{ K) ) 


2 


35.4 


-0.5236 


-6.35 x 10 _a 


-0.225 


4 


1196.8 


0.2024 


-5.92 x 10~ 5 


-0.071 



It is possible to construct eigenfunctions which give 
negligibly small {T 2 ^)j±f, but such functions are not con- 
sistent with the experimental observations. The eigen- 
functions of the ground doublet can be approximated by 

i>a = C 5 ||> +Cl||> +C 3 | ~ |) , (2) 

1>b = C3||> + ci| - |> + cb| - |) , (3) 

where c\ = cos#, C3 = sin 9 cos <p, and C5 = sin (9 sin tp^S 

These functions diagonalize 2 for even K , and linear 

combinations of (ip a + tpb) / V% and {ip a — tpb)/V^ also di- 

(2) 

agonalize J x . 9 — 7r/4 and ip ~ 0.227T give (0 2 ) and 

(0 2 4 ^) around their maximum values, which are close to 
the values in Table [I] Almost identical charge distri- 
bution with that of Fig. ^| is obtained. The ratio of 
(T 2 (4) ) 4 //(T 2 (2) ) 4 / decreases with increasing if and van- 
ishes around ip = 0.737T. However, in this region of (p, 
(2) 

(0 2 ) itself also becomes very small and the anisotropy 
of the dipolar moment within the ab plane disappears; 
(Jx) and {Jy) become to have almost the same values 
although (J z 2 ) is still small. The charge distribution of 
this eigenfunction becomes cylindrical which is elongated 
along the z axis. 

We consider that the dipolar moment is confined along 
(2) 

the x axis by the 2 type quadrupolar moment because 
the magnetic structure in the AFM phase strongly re- 
flects the underlying AFQ order. Therefore, the region 
around tp = 0.737T is not probable. The wavefunction of 
the ground state should have a rank 4 moment that is 
not negligible compared with the rank 2 moment. The 
charge distribution should be like the one in Fig. 1101 

However, the azimuthal-angle dependence of the RXS 
at the E2 resonance indicate that (T 1 2 4 ' > ) 4 / is much 

smaller than (T 2 ) 4 /. This disagreement could be due 
to the r.m.e in Tabic [I] which is the value calculated for 
the atomic wavefunctions within the framework of the 
idealized scattering length^ Although we consider that 
the rank 4 moment should be reflected in the RXS of the 
E2 resonance, the actual intensity should be re-examined 
by more realistic theoretical calculations which take into 
account complex intermediate states more properly. 
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B. d-f Coulomb and quadrupole-strain interactions 

Here, we simply assume that an electron in Dy with 
a wavefunction extended toward B and C atoms, which 
are negatively charged, costs high energy because of the 
Coulomb interaction. Then, the analysis of the interfer- 
ence effect indicates that the charge distribution of the 
5d electron is determined by the d-f Coulomb interac- 
tion, and not by the lattice distortion of B and C. The 
lattice distortion is caused by the direct interaction with 
the 4/ electrons, not with the 5d electrons although they 
are more extended. This conclusion is contrary to the 
theoretical study of Rcf. El where it is concluded that 
the lattice distortion directly modulates the 5d states and 
the .El resonance does not reflect the 4/ quadrupolar mo- 
ment. Later on, however, the same authors reconsidered 
the mechanism and reported that the El resonance can 
be caused by the d-f Coulomb interaction^ 

In the course of the analysis, we assumed that F\ at and 
the El and E2 resonances interfere coherently; the mix- 
ing parameters r± and r-i were fixed at unity, a positive 
real number. Here, it should be noted that we approxi- 
mate the energy spectrum of the resonance by only two 
oscillators at A^i and Ae2, which is actually composed 
of many oscillators of the respective intermediate states. 
In this sense, the mixing parameter should be in gen- 
eral a complex number. Nevertheless, in view of the fact 
that all the experimental data, including the (0 2.5) re- 
flection examined in the previous paper^^i are explained 
consistently by assuming the mixing parameters as posi- 
tive real numbers, we consider that the above conclusion 
must be the case. This point also should be re-examined 
by a more realistic theory which properly considers the 
intermediate states^ 



V. CONCLUSIONS 

We have performed RXS for the (3 1.5) superlattice 
reflection in DyB2C2. Resonant enhancement was ob- 
served below Tq, indicating the AFQ order in which the 
principal axis of the quadrupolar moment is canted in 
accord with the site symmetry of 2/m. 

Azimuthal-angle dependence of the intensities for the 
El and E2 resonances have been measured precisely. We 
analyzed the results using a theory which directly con- 
nects atomic tensors with the scattering amplitudes of 
RXS. The fitting results indicate that the data are well 
explained by the atomic tensors up to rank 2. The rank 
4 tensor is little reflected in the data. However, a simu- 
lation of the wavefunction in the AFQ phase predicts a 
relatively large value of the rank 4 moment in compar- 
ison with that of rank 2. There should be some reason 
the rank 4 moment is difficult to observe in the E2 pro- 
cess of RXS, which should be examined by more realistic 
theoretical calculations properly taking into account the 
complex intermediate states. 

Finally, an intra-atomic d-f Coulomb interaction and a 



direct coupling between the 4/ quadrupolar moment and 
the lattice have been evidenced by the analysis of the in- 
terference structure among the Thomson scattering from 
the lattice distortion, the E2 resonance which observe 
the 4/ state, and the El resonance which observe the hd 
state. The charge distributions of the 5d and 4/ electrons 
are extended along the directions almost perpendicular 
with each other. The B and C atoms along the direction 
where the 4/ charge distribution is extended move away 
from Dy, while those along the direction where the 5<i 
charge distribution is extended move toward Dy. 
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VII. APPENDIX 

We use a formalism developed by Lovesey et al. to 
analyze the experimental results of RXS. This formalism 
simply deals with the intermediate states as an atomic 
core hole, and then the resultant scattering amplitude is 



directly connected with an atomic tensor (T< 



(*)\ 



The 



atomic scattering factor for a resonant process between 
core and valence states with orbital quantum numbers I 
and I, respectively, is expressed as 

K ( ' 



x£<TW>xW(-l)« 
i 

Z (E2) = B J2(-l) K V2K + l{ J K x \ 



(4) 



(5) 



where A = {l\\C(l)\\l)(l\\C(l)\\l){l\R$) 2 and B 



yf{(l\\C(2)\\l)(l\R 2 \l)} 2 are the constant factors for El 



and E2 processes, respectively. Xq and H q ^' are the 
spherical tensors which express the state of x rays in- 
volved in the scattering process: 



r(K) 



x( q K) = E <,V<WI*?> ( 6 ) 
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and expressed in the form of spherical tensors of rank 1. 
1 

h„ = eukw{luli/\2n) (7) 

u,u f =— 1 

2 The structure factor for the (3 1.5) reflection is cal- 

fff) = ti^{2^'\Kq) , (8) culated by F = Z(l) - Z{2) - Z(3) + Z(4), where the 

u,u'=— 2 index represents the four Dy ions in the unit cell. The 

final scattering amplitude at resonance is proportional to 

where e (e') and k (k') are the polarization and wavevec- F/(E — A + iT/2), where A and T are the energy and 

tors of the incident (scattered) x ray, and they are also width of the resonance. 
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